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ADAPTIVE  PHASE-SHIFTER  NULLING  TECHNIQUES  FOR 
LARGE-APERTURE  PHASED  ARRAYS 


1.0  INTRODUCTION 

The  particular  type  of  antenna  system  addressed  is  a  large  aperture  (in  wavelengths)  phased 
array  of  low-sidelobe  design  where  the  investment  is  already  considerable  and  one  simply  could  not 
afford  to  make  it  fully  adaptive.  The  Rome  Air  Development  Center  (RADC)  space-based  radar 
(SBR)  antenna  system  concept  [1]  is  an  excellent  application  case  in  point.  The  low  sidelobes  are 
necessary  for  good  initial  protection  against  sidelobe  interference  and  large  clutter  returns.  Also,  they 
reduce  the  problems  associated  with  strong  sources,  i.e.,  the  number  of  adaptive  degrees-of-freedom 
(DOF)  required  and  the  dynamic  range  of  the  adaptive  processor.  Thus,  retention  of  the  low 
sidelobes  is  considered  a  major  requirement  for  the  system. 

Several  partially  adaptive  array  techniques  are  suggested  in  the  literature  for  approaching  this 
problem,  including  the  conventional  sidelobe-canceller  (SLC)  configuration  |2 ,3) .  constrained  beam- 
space  adaptive  systems  [4-8],  overlapped  subarray  (space-fed  lens)  configurations  [1.9-11],  and 
phase-shifter  nulling  techniques  112-14],  A  brief  commentary  on  each  of  these  solutions  follows. 

1.1  Conventional  SLC 

Figure  1  illustrates  the  conventional  SLC  configuration  consisting  of  a  high-gain  mainbeam 
antenna  surrounded  by  an  auxiliary  array  of  a  few  low-gain  elements.  This  configuration  has  proven 
popular  over  the  years  because  it  can  provide  modest  protection  against  sidelobe  interference  for  the 
least  cost,  with  the  cost  being  roughly  proportional  to  the  number  of  DOF  implemented.  For  our  pur¬ 
poses,  the  SLC  concept  is  economically  attractive  as  a  starting  point;  however  the  specific  arrange¬ 
ment  of  Fig.  1  has  drawbacks  because  the  auxiliary  array  is  separate  from  the  mainbeam  antenna  and 
because  it  operates  in  e'ement  space.  These  make  it  virtually  impossible  to  achieve  low  adaptive 
sidelobes  because  of  the  inherently  wide  element  pattern  coverage,  the  presence  of  array  configuration 
grating  lobes,  and  the  differing  polarization  responses.  Consider  the  two  SLC  auxiliary  pattern 
options  indicated  in  Fig.  2,  where  we  show  cancellation  of  a  single  emitter  by  means  of  a  single  aux¬ 
iliary  element  (dashed  lines)  and  a  single  auxiliary  beam  (solid  lines).  Note  that  the  emitter  is  nulled 
out  by  either  option,  but  the  adapted  pattern  sidelobe  effects  are  quite  different.  The  auxiliary  SLC 
beam  has  little  effect  on  the  remainder  of  the  adapted  pattern  sidelobes  because  of  its  highly  localized 
interaction  beamwidth  region,  whereas  the  auxiliary  SLC  element  interacts  across  the  entire  adapted 
pattern  sidelobe  region  because  of  its  broad  element  pattern.  In  general,  the  auxiliary  elements  option 
produces  increasing  adaptive  sidelobe  degradation  as  the  number  of  interference  sources  increases. 
Also,  there  can  be  a  significant  reduction  in  output  signal-to-noise  ratio  (SNR)  because  of  coupling-in 
additional  receiver  noise  along  with  the  auxiliary  element  signals  [2J. 

1.2  Constrained  Beamspace  Antenna  Systems 

Many  different  possible  configurations  are  within  this  category.  Figure  3  is  a  schematic  diagram 
of  one  particular  type  of  constrained  beamspace  system  [15].  It  consists  of  two  beamformer  subsys¬ 
tems  that  connect  into  common  array  aperture  elements.  The  use  of  common  elements  helps  to  avoid 
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Fig  2  —  Conceptual  cancellation  comparison  between  a  single 
auxiliary  element  and  a  single  auxiliary  beam 
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Fig.  3  —  Constrained  boamspace  SLC  system  with  a  tapped  delay  line  on  each  assigned  beam 


polarization  problems  in  the  adaptive  nulling  performance.  The  mainbeam  subsystem  (on  the  left- 
hand  side  in  Fig.  3)  is  intended  to  function  as  a  conventional  low-sidelobe  phased  array  capable  of 
electronic  scan  over  the  region  of  interest.  Sidelobe  level  is  determined  by  the  usual  quality  of  the 
elements,  phase-shifter  modules,  and  corporate  feed.  The  second  beamformer  subsystem  is  intended 
to  furnish  a  set  of  SLC  beams  for  selective  subtraction  from  the  mainbeam;  thus  it  is  referred  to  as 
the  SLC  subsystem.  This  SLC  subsystem  is  auxiliary  to  the  mainbeam  subsystem  and  consists  of  a 
beamformer  that  is  coupled  into  the  aperture  at  the  elements  prior  to  the  transmit-receive  (TR)  module 
phase  shifters.  This  beamformer  does  not  require  low-sidelobe  design  and  may  consist  of  either  a 
Butler  matrix  type  or  a  lens  type,  such  as  a  Rotman  lens.  From  the  standpoint  of  potential 
bandwidth,  a  favorable  feature  is  that  all  beams  have  the  same  phase  center,  i.e.,  the  geometric  center 
of  the  array  aperture.  Orthogonality  in  the  family  of  beams  is  desirable  but  need  not  be  precise  for 
our  purposes. 

The  output  beam  ports  connect  into  a  “Beam  Assignment  Selector”  where  they  are  electroni¬ 
cally  switched  into  the  SLC  algorithm  processor  based  on  source  location  estimates.  The  idea  is  to 
connect  only  enough  SLC  beams  to  cancel  a  given  source  distribution  situation,  thus  minimizing 


3 


W  F  GABRIEL 


sidelobe  degradation.  The  general  principles  of  beamspace  SLC  are  described  in  the  literature  [4]  and 
offer  advantages  of  a  stable  mainbeam.  retention  of  low  sidelobes,  very  fast  adaptive  response,  and  no 
adaptive  grating  lobes. 

Each  of  the  selected  beams  feeds  into  a  tapped  delay  line,  usually  referred  to  as  a  transversal 
filter  (TF).  This  arrangement  permits  adaptivity  in  the  time/frequency  domain  [15,16]  in  addition  to 
adaptivity  in  the  spatial  domain.  The  TF  output  taps  then  feed  into  the  black  box  labeled  “SLC  Algo¬ 
rithm  Signal  Processor,’’  which  applies  an  adaptive  algorithm  to  obtain  the  TR  tap  weights  for  achiev¬ 
ing  cancellation.  This  type  of  system  is  amenable  to  any  of  the  current  adaptive  processing  algo¬ 
rithms,  including  analog  versions. 

Figure  4  shows  typical  patterns  obtained  from  computer  simulations  for  a  16  element  linear 
array.  The  mainbeam  subsystem  is  given  a  quiescent  Taylor  illumination  designed  for  30  dB 
sidelobes,  and  a  Butler  matrix  was  chosen  for  the  SLC  subsystem  beamformer.  A  two-beam  cluster 
consisting  of  SLC  beams  numbers  10  and  1 1  were  assigned  to  cover  two  strong  40  dB  noncoherent 
sources  located  at  18°  and  22°  azimuth.  Note  that  the  adapted  pattern  shows  cancellation  of  both 
sources  with  little  perturbation  of  the  mainbeam  sidelobes,  except  in  the  immediate  vicinity  of  the 
sources. 

Figure  3  indicates  that  a  third  subsystem  is  the  source  estimation  processor.  This  is  an  all- 
digital  processor  that  operates  on  the  signals  received  from  a  few  of  the  aperture  elements  to  detect 
and  track  the  interference  sources  of  interest.  Processing  the  digital  signals  to  estimate  the  sources 
may  be  accomplished  with  a  number  of  spectral  estimation  algorithms  available  in  the  literature  [4] 
Once  the  source  estimation  information  is  available,  we  assign  beamformer  beams  by  means  of  a 
computer-logic  program.  It  is  important  to  note  that  beam  assignment  does  not  require  great  accu¬ 
racy,  i.e.,  a  half-beam  width  is  usually  close  enough.  Also,  clusters  of  two  or  three  adjacent  beams 
may  be  assigned  for  doubtful  cases. 

1.3  Overlapped  Subarray  Antenna  Systems 

As  in  the  previous  category,  there  are  many  possible  configurations  within  overlapped  subarray 
antenna  systems.  The  type  chosen  for  illustration  is  taken  from  Southall  and  McGrath  [1],  where  Fig. 
5  shows  the  antenna  layout  drawing  for  a  completely  overlapped  subarray  antenna  (COS A).  It  con¬ 
sists  of  a  space-fed,  two-dimensional  bootlace  lens  system  in  which  the  space-feed  is  a  Rotman  lens 
beamformer  [17].  Figure  6(a)  illustrates  excitation  of  one  of  the  beamformer  inputs,  or  subarray 
ports,  resulting  in  a  sin(x)/x  illumination  of  the  bootlace  lens.  The  bootlace  lens,  in  turn,  then 
transforms  its  illumination  into  a  near-rectangular  far-field  pattern  function  that  may  be  scanned  by 
phase  shifters  behind  the  objective  len's  front  face.  Figure  6(b)  illustrates  how  all  of  the  subarray 
ports  may  be  driven  simultaneously  (the  resulting  illuminations  on  the  objective  lens  completely  over¬ 
lap  one  another)  and  are  coherently  combined  and  weighted  to  synthesize  a  tapered  illumination  on  the 
objective  lens,  thereby  producing  a  low  sidelobe  far-field  pattern.  Because  of  the  Fourier  transform 
principles  of  operation,  this  type  of  space-feed  is  often  referred  to  as  a  transform  feed. 

A  distinct  advantage  permitted  by  this  type  of  antenna  system  is  wide  instanteous  bandw  idth  [9], 
obtained  by  applying  the  appropriate  time  delays  at  subarray  ports.  Note  that  the  number  of  time 
delays  required  in  the  transform  feed  is  far  less  than  the  number  that  would  be  required  to  time-delay 
steer  the  objective  lens,  and  yet  it  can  synthesize  an  approximate  equivalent  time-delay  steering  over 
considerable  bandwidth. 

Since  the  subarray  phase  centers  are  located  many  wavelengths  apart,  the  objective-lens  aperture 
array  factor  contains  grating  lobes  as  illustrated  in  Fig.  7.  The  resultant  overall  far-field  pattern  is  the 
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Fig.  5  —  Antenna  layout  drawing  of  COSA  system  (Ref.  I.  p.  468) 
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SYNTHESIZED 


(a)  (b) 


Fig.  6  —  COSA  fundamental  operating  principles;  (a)  single  subarray  beam, 
and  (b)  weighted  sum  of  subarray  beams  (Ref.  1,  p.  466) 
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Fig.  7  —  COSA  pattern  characteristics;  (a)  /  =  /„; 
(b)  /  >  /0;  and  (c)  /  <  /0  (Ref.  I,  p.  466) 
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product  of  this  array  factor  and  the  flat-topped  subarray  pattern.  The  subarray  pattern  can  be 
considered  as  a  window  that  selects  the  array  factor  mainlobe  at  a  given  steering  direction  and 
suppresses  all  grating  lobes.  System  bandwidth  is  eventually  limited  by  the  fact  that  the  subarray  pat¬ 
tern  is  phase-steered  and  shifts  in  angle  vs  frequency  as  illustrated  in  Fig.  7.  Also,  although  the  array 
factor  mainlobe  remains  fixed  because  of  its  time-delay  steering,  all  grating  lobes  shift  in  angle  vs  fre¬ 
quency  and  may  abet  intrusion  into  the  shifting  subarray  window. 

From  a  partially  adaptive  antenna  point  of  view,  this  configuration  offers  three  basic  adaptive 
control  alternatives: 

•  The  subarray  input-beam  ports  on  the  transform  feed. 

•  The  focal-plane-element  radiators  on  the  transform  feed. 

•  The  objective  lens  element  phase-shifters  and/or  element  transceiver  modules. 

It  is  beyond  the  scope  of  this  report  to  consider  these  alternatives  in  depth  because  of  the  complexity 
of  the  systems.  However,  one  can  make  some  general  observations  based  on  Figs.  5,  6.  and  7.  The 
first  alternative  would  require  the  fewest  adaptive-weight  DOF  and  should  enable  excellent  adaptive 
performance  within  the  confines  of  the  flat-topped  subarray  pattern,  i.e.,  in  the  “near-in”  sidelobe 
region.  Outside  of  the  subarray  pattern  (in  the  “far-out”  sidelobe  region),  adaptive  performance  is 
inferior  because  of  low  gain  conditions  and  a  severe  coupled  disturbance  into  the  “near-in"  sidelobe 
region,  i.e.,  the  entire  subarray  pattern  must  fluctuate  in  accordance  with  adaptive  weight  fluctuations. 
Mainbeam  and/or  sidelobe  level  adaptive  constraints  are  required  to  prevent  excessive  perturbations 
from  "far-out”  sidelobe  interference. 

The  second  alternative  has  performance  characteristics  similar  to  the  first,  since  they  are  related 
by  a  linear  transformation.  However,  some  differences  are  possible  because  of  the  change  in  objec¬ 
tive  lens  illumination  per  adaptive  weight,  i.e.,  each  adaptive  weight  now  controls  a  pencil  beam  in 
the  far  field.  One  could  consider  certain  combinations  of  amplitude-only  and  phase-only  adaptive 
control  that  may  be  more  economical  than  the  first  alternative. 

The  third  alternative  should  enable  the  best  “far-out”  sidelobe  region  adaptive  performance 
because  its  adaptive  DOF  can  be  designed  to  cover  the  entire  sidelobe  region  of  interest,  even  when 
the  objective  lens  aperture  is  divided  into  a  reasonable  number  of  subarray  “elements.”  Mainbeam 
constraints  should  not  be  required.  Note,  however,  that  this  alternative  is  susceptible  to  bandwidth 
degradation  and  the  hardware  cost  will  be  higher  because  of  signal  acquisition  transfer  complexity. 
Signals  from  the  adaptively  controlled  subarray  “elements”  must  be  acquired  and  transferred  from 
the  lens  to  the  receiver-processor.  Then  the  resultant  control  signals  from  the  receiver-processor  are 
sent  back  to  the  lens  and  distributed  among  the  elements. 

1.4  Phase-Shifter  Nulling  Techniques 

The  particular  configuration  within  this  category  chosen  for  discussion  is  patterned  after  Baird 
and  Rassweiler  |I2],  When  combined  with  simple  adaptive  search  algorithms,  this  nulling  method  is 
the  most  simple  and  least  expensive  of  all  adaptive  techniques  available.  It  does  not  require  auxiliary 
elements,  correlators,  or  beamformers  but  depends  on  access  to  the  array  element  phase  shifters  for 
injection  of  small  phase  (and/or  amplitude)  perturbations  of  the  mainbeam-aperture  distribution.  Fig¬ 
ure  8  shows  a  typical  element-vector  perturbation  diagram  for  a  phased-array  antenna  that  uses  ferrite 
phase  shifters.  One  may  view  the  perturbation  as  resulting  from  superposition  of  the  wavefront  for 
the  mainbeam  plus  a  wavefront  for  the  emitter  direction.  Basic  principles  of  aperture  ripple- 


8 


NRl.  REPORT  9132 


modulation  effects  are  described  in  Ref.  18,  and  beamspacc  decomposition  for  low-sideiobe  phased 
arrays  is  described  in  Ref.  12. 

The  small  phase  (and/or  amplitude)  perturbations  illustrated  in  Fig.  8  may  be  injected  by  means 
ot  several  techniques,  depending  upon  the  particular  type  ot  components  and  control 
hardware/software  available: 

•  The  phase  shifter  itself  may  permit  small  analog  vernier  control,  i.e.,  analog  ferrite  phase 
shifters,  or  even  digital  phase  shifters  in  which  one  of  the  bits  may  be  perturbed  by  vernier 
analog  current  drive. 

•  If  the  array  contains  hundreds  or  thousands  of  elements,  the  technique  of  “row/column" 
averaging  may  be  used  to  get  small,  net,  phase  perturbations  in  azimuth  or  elevation,  even 
though  the  individual  phase  shifters  have  only  four-  or  five-bit  control. 

•  The  addition  of  a  vernier  attenuator  or  phase  shifter  at  each  element. 

•  The  existence  of  a  separate,  controlled  beamformer  coupled  into  the  aperture  elements. 


The  adaptive  algorithms  needed  for  adjusting  the  element  weight  perturbations  are  generally 
known  as  'search  techniques”  (13,161,  and  they  involve  using  the  value  of  some  performance  criteria 
(usually  mainbeam  output  power)  to  determine  the  direction  of  the  next  perturbation.  This  class  of 
algorithm  is  usually  characterized  by  relatively  slow  convergence  speed  since  individual  element 
direction  gradients  are  not  available  and  must  be  approximated  by  means  of  time-consuming  random 
trials.  However,  in  the  current  application  where  we  know  the  direction  of  the  emitters,  their  relative 
power  levels,  and  the  mainbeam  characteristics,  beamspace  decomposition  search  techniques  can  steer 


9 


W.  F  GABRIEL 


sidelobe  nulls  onto  the  emitter  positions  with  surprising  conveigence  speed.  This  report  examine:;  this 
latter  technique  in  some  detail  to  gain  a  better  perspective  of  its  performance/cost  tradeoff  characteris¬ 
tics  when  applied  to  large  phased-array  antenna  systems. 

Section  2  develops  the  aperture  ripple  modulation  algorithms  that  were  used  in  this  investigation; 
Section  3  discusses  simulation  results  based  on  those  algorithms;  Section  4  presents  beamspacc 
options  for  improving  performance  at  the  cost  of  the  additional  DOF,  and  Section  5  presents  the  con¬ 
clusions  reached  in  this  study. 

2.0  APERTURE  RIPPLE  MODULATION  ALGORITHMS 

The  beamspace  decomposition  discussed  by  Baird  and  Rassweiler  [12]  assumes  that  the  source 
direction  angles  are  known,  whereupon  we  are  working  with  the  simple  element  vector  diagram  illus¬ 
trated  in  Fig.  8.  This  diagram  represents  the  superposition  of  the  element  phasing  for  our  quiescent 
mainbeam  wavefront,  plus  the  phasing  for  each  of  /  emitter  small-amplitude  wavefronts.  Thus,  by 
referring  to  the  vector  diagram,  a  mathematical  expression  may  be  written  for  the  Ath  e'  ment  weight 
B'(A)  of  a  linear  array. 


mk)  =  wqo<) 


! 

1  +  A,  exp  (i p,(k )) 
i  =  i 


( 1 ) 


where 

Wq(k)  is  the  quiescent  beam  A; th  element  weight. 

A,  is  the  amplitude  of  the  /  th  source  perturbation,  and 

i/-,  (k )  is  the  phase  of  the  /th  source  perturbation  at  the  A  th  element. 

Further.  \p,(k )  is  referenced  to  the  quiescent-mainbeam  direction  element  phasing  and  may  be 
expanded  in  the  form 


'Pi(k)  = 


2  7T(V 

A 


(sin  0, 


-  sin  0O) 


k 


r  ■> 

K  +  1 

l  2  J 

(2) 


where 

d  is  the  element  spacing, 

X  is  the  wavelength, 

K  is  the  total  number  of  elements, 

0O  is  the  quiescent  mainbeam  steering  direction,  and 

0,  is  the  direction  angle  of  the  /  th  source. 

The  array  elements  are  assumed  to  be  equally  spaced,  and  all  element  phasing  is  referenced  to  the 
midpoint,  or  center,  of  the  array. 

The  only  unknown  quantity  in  Eqs.  (I)  and  (2)  is  the  coefficient  A,  of  the  small-amplitude 
source  wavefronts.  These  coefficients  are  arrived  at  by  a  recursive  modulation  trial-and-update  rou¬ 
tine  that  may  be  expressed  in  the  form. 


A,(n  +  1)  = 


Tin ) 

r(n)  +  1 


A,(n) 


+  fijin), 


(3) 
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where 


A,  (n ) 

A, (n  4  1) 
H,  (n ) 
r(n ) 


is  the  valuf  of  A ,  used  in  current  trial, 
is  the  new  value  of  A,  for  the  next  trial, 
is  the  update  increment  from  current  trial,  and 
is  the  integration  time  constant. 


Note  that  the  A ,  are  quantities  that  are  formed  by  integrating  the  update  increments,  and  that  they 
incorporate  a  decay  factor  determined  by  the  time  constant  r(n).  This  time  constant  is  chosen  to  have 
a  SNR  dependency  as  follows. 


T(n )  =  7.5  + 


80 
(« ) 


(4) 


where  S(n)  is  the  current  trial  value  of  best  SNR. 


The  update  increments  are  determined  from  a  weight  perturbation  modulation  trial  that  involves  a  sim¬ 
ple  plus-minus  sequence,  with  the  array  output  power  P(t),  monitored  and  averaged  for  equal  time 
periods  during  the  modulation  trial.  The  sequence  proceeds  as  follows: 

Step  1:  If  this  is  a  new  start,  a  value  of  averaged  P(i)  for  unperturbed  weights  W  is  required. 

Let  P0  represent  that  measured  average  power  level.  Note  that  W  may  be  the  quies¬ 
cent  mainbeam  weights  if  the  initial  A,  are  set  to  zero  for  a  new  start  (Eq.  (1)). 

Step  2:  If  this  is  a  new  start,  an  initial  value  for  the  weight  perturbation  magnitudes  A,  is 
required.  Initial  magnitudes  are  usually  chosen  about  6  dB  below  the  quiescent  beam 
sidelobe  level.  For  example,  if  the  nominal  sidelobe  level  is  -30  dB,  then  we  would 
set  the  initial  values  of  all  A,  =  0.017,  or  -36  dB  level. 

Step  3.  The  A,  are  applied  as  a  weight  perturbation  modulation  in  plus-minus  sequence  for  the 
/th  source  alone,  where 

W,  +  (k)  =  IF(*)[1  +  A,  exp  (*/(*))],  (5) 


and 

IT,-!*)  =  -  A,  exp  (\pj (*))],  (6) 

where  was  given  in  Eq.  (2).  The  output  power  is  monitored  for  this  sequence, 
giving  us  two  averaged  power  levels, 

P,{  +  )  is  the  averaged  P{t)  for  W,+  weights,  and 
P,(~)  is  the  averaged  P(t)  for  W,"  weights. 


Note  that  the  phase  term  \pj(k)  confines  the  modulation  to  a  particular  source;  we 
are  essentially  adding  and  subtracting  a  low-level  beam  pointed  in  the  direction  of  the 
/  th  source  alone.  This  is  a  heamspace  low-level  modulation  technique. 
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Step  4:  Determine  best  A,  from  the  modulation  trial  of  Step  3.  The  following  logic  procedure 
is  used  in  a  Fortran  IV  software  program: 

F  =  0 

IF  ((P,(  +  )  >  P0).  AND.  (/>  (-)  >  p0))  F  =  1  ;  set  flag. 

IF  (F  =  1)  GO  TO  33  ;  source  in  null,  do  not  change. 

u  =  +1 

X  =  Pj(  +  )  —  Pj(-)  ;  modulation  difference  in  power  level. 

IF  ( X  >  0)  GO  TO  10 

Z  =  Po/Pj(  +  )  ;  power  level  change  ratio. 

n  n  A 

'  o  =  ‘  /(+)  ;  Po  ratchet-down  to  lowest  power  level. 


GO  TO  15 


10  u  =  -1 


Z  =  p0/p,(~) 

Pq  ~  Pi(~ ) 

15  Hj(rt)  =  u  *  A, 

S(n)  =  Pq/Nq 

r(n)  =  7.5  +  (80 /'^S(n)) 


;  power  level  change  ratio. 

;  P0  ratchet-down  to  lowest  power  level. 
;  retain  the  best  A, . 

;  current  best  SNR. 

;  new  time  constant. 


With  a  new  value  for  n,(n)  and  r(n),  we  can  update  coefficient  A,  (n)  per  Eq.  (3), 
and  therefore  W(k)%  per  Eq.  (1). 

Step  5:  The  magnitudes  of  the  A,  do  not  remain  constant.  They  must  diminish  as  the  output 
power  level  P0  is  driven  toward  the  receiver  noise  level;  otherwise  the  modulation 
perturbations  will  overdrive  the  desired  null  conditions  and  result  in  poor  convergence 
behavior.  The  following  logic  procedure  is  used  to  continue  from  Step  4. 
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X  =  7.5/-(rt) 

Y  =  10  *  ALOGIO  (Z) 
B  =  10 /(Y  +  10) 

Q  =  0.022  *  X  *  B 
DO  20  J  =  1 ,  / 

IF  ((2  >  A,  )  GO  TO  20 

A,  =  0 


;  dependence  on  current  best  SNR. 
;  power  level  in  dB. 

;  dependence  on  power  level  drop. 
;  maximum  allowable  value. 


;  reset  A,  to  new  lower  value. 


20  CONTINUE 

Note  that  since  P0  depends  on  all  source  nulling,  then  all  A,  are  set  equal  except 
for  retention  of  a  smaller  prior  value.  This  permits  a  ratchet-down  convergence 
behavior. 

Step  6:  If  a  given  source  is  already  in  a  sidelobe  null,  then  the  flag  (F  =  1)  set  in  Step  4  is 
used  to  bypass  any  change  in  the  A,  coefficients  of  Eqs.  (1)  and  (3).  However,  the 
modulation  magnitude  A,  for  that  particular  source  is  reduced  recursively  to  a  very 
small  value,  essentially  deleting  it  from  the  processing.  The  Fortran  IV  software 
logic  statement  entered  from  Step  4  is  as  follows: 

33  IF  (F  =  1)  A,  =  0.5  *  A,  +  0.00035 

In  this  statement,  the  value  of  A,  is  reduced  by  6  dB  on  each  flag  (F  =  1)  recursion, 
reaching  an  ultimate  minimum  of  —63  dB,  if  the  source  continues  to  remain  in  a 
sidelobe  null. 

Step  7:  If  a  given  emitter  happens  to  turn  off.  this  would  be  detected  by  the  source 
tracking/estimation  subsystem,  and  it  would  soon  be  deleted  from  the  high-priority 
source  directory  kept  in  memory.  Thus,  no  further  beam  modulations  would  be  made 
in  its  direction. 

2.1  Modifications  for  Phase-Only  Nulling 

The  algorithm  described  incorporates  both  amplitude  and  phase  in  the  weight  perturbations  and 
will  result  in  the  best  overall  nulling  performance.  However,  most  phased  array  applications  would 
require  restriction  of  the  perturbations  to  phase-only  changes,  so  we  introduce  the  modifications 
described  here. 

First,  with  reference  to  the  vector  diagram  of  Fig.  8,  only  the  phase  angle  change  is  applicable 
so  that  Eq.  (1)  becomes, 
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W(k)  =  Wq(k) 


i  +j  E  Ai sin  (iM*)) 

i  —  1 


(7) 


where  it  is  noted  that  only  the  imaginary  part  of  the  A;  terms  is  utilized.  This  permits  us  to  under¬ 
stand  the  interesting  phenomenon  of  unwanted  “symmetrical  lobes”  that  appear  in  phase-only  modu¬ 
lation,  because  we  can  substitute  the  following  trigonometric  identity  into  Eq.  (7), 

j  sin  (*,-(*))  =  y  -  e~j*M)\.  (8) 

Obviously,  we  have  two  lobes  located  at  ±ipj(k),  and  from  Eq.  (2)  we  see  that  they  are  symmetrica) 
with  respect  to  the  pointing  direction  of  the  mainbeam.  An  example  pattern  plot  of  such  lobes  is  dis¬ 
cussed  in  Section  3. 

Equations  (2)  through  (4)  remain  unchanged,  but  in  Step  2  it  is  necessary  to  double  the  initial 
value  for  the  weight  perturbation  magnitudes  A,  .  For  the  example  given,  we  now  set  all  A,  =  0.034, 
or  —30  dB  level.  The  reason  for  doubling  is  that  phase-only  modulation  is  roughly  half  as  effective 
as  the  combined  amplitude  plus  phase. 

In  Step  3,  we  again  encounter  the  phase  angle  restriction  so  that  Eqs.  (5)  and  (6)  become 

W^ik)  =  W(k)  [1  +  j  A,  sin  (i £,■  (*))],  and  (9) 

W~(k)  =  W(k)  11  -  j  A,  sin  (*,•(*))].  (10) 

Steps  4  and  5  remain  the  same,  but  in  Step  6  we  again  double  the  minimum  modulation  value  for  the 
flag  (F  =  1)  condition, 

33  IF  (F  =  1)  A,  =  0.5  *  A,  +  0.0007 


such  that  the  ultimate  minimum  value  is  now  -57  dB  if  the  source  continues  to  remain  in  a  sidelobe 
null. 

2.2  Other  Algorithms 

The  algorithm  just  described  is  a  beamspace,  single-cycle,  plus-minus  perturbation  sequence 
technique  that  uses  the  changes  in  smoothed/averaged  output  power.  There  are  also  other  algorithms 
and  techniques  described  in  the  literature  that  would  be  suitable  for  our  application.  For  example,  the 
orthogonal  binary  code  sequences  described  by  Cantoni  and  Godara  1 19,20]  could  be  adapted  to  our 
situation,  where  it  is  noted  that  they  correlate  the  instantaneous  output  power  with  their  sequences 
rather  than  the  smoothed  output  power.  Also,  the  coherent  perturbation-algorithm  work  of  Davis  et 
al.  [21]  should  be  applicable.  These  other  techniques  are  not  pursued  in  this  investigation. 

3.0  SIMULATION  RESULTS  FROM  APERTURE-RIPPLE  ALGORITHM 

In  this  section,  we  describe  performance  results  obtained  by  applying  the  algorithms  of  Section  2 
to  computer-simulated  source  situations.  Figure  9(a)  illustrates  an  example  wherein  we  have  two 
rather  strong  emitters  of  46  dB  SNR  located  in  the  sidelobe  region  of  a  16  element  linear  array  with 
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Fig.  9  —  Quiescent  mainbeam  and  aperture  ripple-modulation  beams  for  16  element 
linear  array,  two  46  dB  sources  at  10°  and  42°;  (a)  quiescent  mainbeam  pattern.  30 
dB  Taylor  weighting,  and  (b)  subtraction  beams  resulting  from  phase-only 
perturbation 


half-wavelength  element  spacing,  and  a  30  dB  Taylor  illumination.  The  sources  are  noncoherent  and 
are  swept  in  frequency  over  a  5%  RF  bandwidth.  Since  we  know  the  directions  of  the  sources,  phase 
perturbations  are  applied  across  the  array  to  generate  two  beams  in  those  directions.  Figure  9(b) 
illustrates  pattern  plots  of  the  two  beams,  labelled  1  and  2;  note  that  these  beams  are  the  result  of 
phase-only  perturbation.  They  have  good  sidelobe  levels  because  they  automatically  incorporate  the 
30  dB  taylor  illumination  taper,  but  they  are  also  subject  to  the  appearance  of  unwanted  “symmetrical 
lobes”  because  of  the  phase-only  approximation  discussed  in  Section  2. 
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The  symmetrical  lobe  appears  on  the  opposite  side  of  the  mainbeam,  i.e.,  the  mainbeam  pointing 
direction  is  the  point  of  symmetry  as  seen  in  Fig.  9  for  the  Number  1  beam.  Note  in  Fig.  9  that  the 
symmetrical  lobe  for  the  Number  2  beam  does  not  appear  in  visible  space  because  it  is  too  far  beyond 
the  mainbeam  direction.  If  the  perturbations  were  to  incorporate  both  amplitude  and  phase,  then  the 
unwanted  symmetrical  lobes  would  not  be  formed. 

As  discussed  in  section  2,  the  search  algorithm  applies  a  starting  perturbation  magnitude  level 
that  is  roughly  equal  to  the  mainbeam  sidelobe  level,  and  it  monitors  the  averaged  power  output  for  a 
plus-minus  sequence  modulation.  If  the  power  output  drops  below  the  P0  level  as  measured  on  the 
previous  trial,  that  perturbation  is  retained  and  the  value  of  P0  is  reset  to  the  new  lower  level.  Figure 
10(a)  shows  the  resultant  “ratchet-down”  effect  on  the  output  power  level,  where  it  is  noted  that  con¬ 
vergence  occurs  irt  about  600  "snapshots”  for  this  case. 

A  snapshot  is  equivalent  to  the  element  data  contained  in  one  range  bin  sampled  simultaneously 
across  the  array  [4|.  In  Fig.  10.  the  dot  symbol  ".”  denotes  the  instantaneous  power  output  of  the 
mainbeam  from  each  snapshot;  the  clusters  of  “  +  "  symbols  denote  an  average  of  the  instantaneous 
power  taken  over  the  previous  21  snapshots;  they  clearly  identify  the  P,l  +  )  and  P'  ~ '  power  levels 

associated  with  a  given  weight  perturbation.  The  horizontal  line  symbols  " _ ”  denote  the 

unperturbed  average  power  level  P0  from  the  previous  plus-minus  perturbation  sequence;  the  heavy, 
dashed-line  has  been  added  to  allow  the  reader  to  follow  the  approximate  convergence  of  the  average 
output  power  level  down  to  receiver  noise  level  (0  dB). 

The  reader  may  find  it  helpful  to  go  through  the  first  complete  plus-minus  sequence  indicated  in 
Fig.  10(a).  Starting  from  the  left,  the  initial  new  start  P0  unperturbed  level  is  about  29.5  dB.  The 
first  plus-minus  cycle  is  then  directed  toward  source  Number  1  located  at  + 10°  azimuth;  the  resultant 
power  levels  are  identified  as  P,(  +  )  at  28.5  dB  and  P,(-)  at  33.5  dB.  from  which  the  (  +  )  perturba¬ 
tion  is  retained  and  PQ  is  reset  to  equal  P,(  +  )  at  28.5  dB.  The  second  plus-minus  cycle  is  men 
directed  toward  source  Number  2  located  at  +42°  azimuth;  the  resultant  power  levels  are  identified 
as  P2(  +  )  at  26  dB  and  P 2(-)  at  29  dB,  from  which  the  (  +  )  perturbation  is  retained  and  PQ  is  reset 
to  equal  P2(  +  )  at  26  dB.  The  sequence  is  then  repeated  until  the  output  converges  close  to  the 
receiver  noise  power  level.  Figure  10(b)  illustrates  the  somewhat  faster  convergence  obtained  when 
the  perturbations  incorporate  both  amplitude  and  phase. 

Figure  1 1  shows  a  typical  adapted  pattern  after  convergence.  The  adapted  pattern  remains  quite 
stable  and  has  successfully  nulled  out  both  sources.  A  modest  increase  in  a  sidelobe  can  be  seen  on 
the  left-hand  side  of  the  mainbeam  at  about  -44°,  which  is  the  exact  position  of  the  unwanted 
symmetrical  lobe  discussed  in  Fig.  9(b). 

4.0  BEAMSPACE  TECHNIQUES  FOR  IMPROVING  PERFORMANCE 

Thus  far,  we  have  considered  only  the  most  simple  configuration  where  there  is  just  one  main- 
beam  output  port  to  work  with.  In  this  Section  we  discuss  the  concept  of  having  available  additional 
output-beam  ports.  For  example,  let  us  assume  that  our  phased-array  antenna  system  has  monopulse 
tracking  capability  with  both  a  sum  (£)  and  a  difference  (A)  output  port,  as  shown  in  Fig.  12.  This 
constitutes  an  advantage  for  our  simple  adaptive  system  because  the  availability  of  any  extra  beam 
port  essentially  gives  us  an  auxiliary  beam,  an  extra  DOF.  Because  of  this,  several  different  ways 
exist  for  improving  our  adaptive  performance. 

First,  one  could  simply  apply  the  aperture  ripple  modulation  technique  of  Section  2  separately  to 
the  £  and  A  beams.  This  would  require  superposition  of  the  two  modulations  of  each  element 
weight,  but  the  fact  that  the  £  and  A  beams  are  orthogonal  (or  nearly  so)  should  permit  the  technique 
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SPATIAL  ANGLE  IN  DEGREES 

Fig.  1 1  —  Adapted  pattern  after  convergence  for  the  two-source  case 
of  Fig.  9:  five  independent  trials  plotted 


Y0M  OUTPUT 


Fig  12  —  Typical  monopulse  sum  and  difference  beamformer  network,  with 
adaptive  weight  option  added  for  using  the  difference  beam  as  an  auxiliary 
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to  develop  overlapping  nulls  in  the  two  patterns.  Alternatively,  we  could  work  with  the  two  halves  of 
the  array  before  the  signals  are  combined  into  the  £  and  A  channels.  Haupt  [22]  discusses  the  ramifi¬ 
cations  of  simultaneous  nulling  in  £  and  A  channels. 

Second,  by  reserving  the  A  beam  for  use  as  an  auxiliary  canceller  beam,  one  can  separate  the 
steering  of  a  low-level  lobe  toward  an  emitter  from  the  adaptive  weight  needed  to  subtract  it  from  the 
£  beam.  This  option  is  indicated  by  the  dashed-line  configuration  in  Fig.  12,  and  it  would  have  the 
advantage  of  very  fast  cancellation  for  a  single  interference  source.  The  disadvantages  to  the  option 
include  rather  large  changes  in  the  resultant  adapted  £  pattern;  not  only  are  there  sidelobe  perturba¬ 
tions  in  the  vicinity  of  the  emitter,  but  the  mainlobe  suffers  distortion  because  of  the  subtraction  of  the 
A  beam  from  it.  Note  that  this  option  results  in  the  interesting  phenomenon  that  a  “£  type"  of  ripple 
beam  exists  in  the  A  channel,  and  a  "A  type"  of  ripple  beam  exists  in  the  £  channel. 

Next,  let  us  assume  that  our  phased-array  antenna  system  consists  of  a  space-fed  lens  in  which 
the  mainbeam  is  scanned  by  phase-shifters  located  in  the  objective  lens,  and  we  wish  to  consider  the 
option  of  using  an  o'  f-axis,  defocused  feed  as  an  auxiliary  beam.  Figure  13  illustrates  this  configura¬ 
tion.  The  possibility  here  is  that  aperture  ripple  modulation  could  be  used  to  focus  the  defocused  feed 
onto  the  emitter,  and  then  adaptively  weight  its  output  to  cancel  the  interference  from  the  mainbeam. 
Again,  this  configuration  would  have  the  advantage  of  separating  the  steering  of  a  low-level  lobe  from 
its  adaptive  weight  and  could  achieve  very  fast  cancellation  of  the  interference  source.  Some 
deterioration  in  the  resultant  adapted  pattern  sidelobes  should  be  expected  because  of  the  double 
defocusing  action,  i.e.,  the  mainbeam  sidelobes  will  be  perturbed  by  the  defocused  aperture  ripple, 
and  the  defocused  feed  will  have  a  "major  lobe"  caused  by  “looking  thru"  the  mainbeam  focusing. 
This  major  lobe  appears  in  the  sidelobes  of  the  resultant  adapted  pattern.  Note  that  linear  superposi¬ 
tion  allows  us  to  add  several  defocused  feeds  to  this  type  of  phased-array  system,  thus  making  it 
amenable  to  very  fast  adaptive  response  with  several  DOF. 


APERTURE 

APERTURE  ARRAY  RIPPLE 


Fig  13  —  Space-fed  lens  configuration  utilizing  off-axis,  defocused  feed 
to  obtain  an  auxiliary  beam  for  sidelobe  nulling 
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5.0  CONCLUSIONS 

Several  simple  phase-shifter  nulling  techniques  have  been  investigated  for  application  to  large- 
aperture  phased  arrays.  The  intended  purpose  is  to  provide  a  modest  partial  adaptive  capability" at 
reasonable  cost.  These  techniques  depend  on  having  access  to  the  array-element  phase  shifters  for  the 
injection  of  small  phase  (and/or  amplitude)  perturbations  of  the  mainbeam  aperture  distribution  but 
they  do  not  require  auxiliary  elements,  correlators,  or  beam-formers.  A  key  factor  to  their  successful 
operation  is  a  priori  knowledge  of  the  directions/locations  of  emitter  sources  to  be  nulled  out,  and  this 

must  be  provided  for  either  by  mainbeam  tracking  of  the  sources  or  by  a  separate  source  estimation 
processor  subsystem. 

A  complete  description  is  given  for  the  particular  aperture  ripple  modulation  algorithms  used, 
based  on  beamspace  decomposition  and  adaptive  “search”  techniques.  These  algorithms  are  then 
applied  to  a  computer-simulated  source  situation  involving  two  strong  noncoherent  emitters  that  are 
swept  in  frequency  over  a  5%  bandwidth.  Both  sources  are  nulled  out  within  about  600  snapshots  of 
data,  which  is  a  remarkably  fast  convergence  for  such  a  simple  phase-only  adaptive  system.  The  per¬ 
formance  results  achieved  indicate  that  the  concept  is  feasible  and  should  be  investigated  further  in  an 
experimental  phased-array  system.  Some  additional  beamspace  options  are  also  described,  which 
would  permit  improvements  in  adaptive  performance  if  additional  DOF  were  available. 
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